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Asymptotic stability of wave equations coupled by 

velocities 


Yan Cui* Zhiqiang Wang^ 
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Abstract 

This paper is devoted to study the asymptotic stability of wave equations with con¬ 
stant coefficients coupled by velocities. By using Riesz basis approach, multiplier method 
and frequency domain approach respectively, we find the sufficient and necessary condi¬ 
tion, that the coefficients satisfy, leading to the exponential stability of the system. In 
addition, we give the optimal decay rate in one dimensional case. 
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1 Introduction and Main Results 

In this paper, we consider the long time behavior of the solution to a system of wave 
equations with constant coefficients coupled by velocities. In particular, we want to study 
what kind of conditions, that the coefficients satisfy, lead to the exponential stability of the 
system. 

In the case of scalar wave equation, there are numerous results on asymptotic stability 
or stabilization with internal or boundary damping. Cox and Zuazua [7] studied, by Fourier 
analysis, the energy decay of 

/ 

utt - Uxx ^ a{x)ut = ^ in (0,1) X (0, -boo), 

< u{0,t) = u{l,t) =0 in (0,-|-oo), 

(u,ut)(0) = in (0,1) 
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with indefinite damping. As a corollary, the exponential decay of the energy and its optimal 
rate were shown in [7] when a > 0 in (0,1). Using Multiplier method, Alabau-Boussouira 
Ell! and Alabau et al. [2] proved the indirect stabilization of wave systems coupled by 
displacements, for instance, 

utt — Au + a{x)ut + h{x)v = 0 in U x (0, +oo), 

vtt — Au — b{x)u = 0 in n X (0, +oo), 

< 

u = V = 0 on r X (0, +oo), 

(M,«t)(0) = (u,ut)(0) = in 0 

\ 

That is, the damping is acted only in one equation and the total energy of the whole system 
decays polynomially due to the coupling effect. Using the criteria of polynomial decay in 
Liu and Rao proved in HD, by frequency domain approach, the polynomial stability of a 
partially damped wave system with weak coupling by displacements and multiple propagation 
speeds. Liu and Rao [T8| also proved the indirect stabilization, by Riesz basis approach, with 
optimal polynomial decay rate for a wave system which is coupled by displacements in one 
dimensional case. 

Recently, Alabau-Boussouira, Wang and Yu [S] obtained the indirect stabilization, com¬ 
bining multiplier method and weighted energy techniques, for the following damped wave 
system with variable coefficients coupled by velocities 

Utt — Am -|- p{x, ut) -|- b{x)vt = 0 in n X (0, -|-oo), 

Vtt — Am — b{x)ut = 0 in n X (0, -|-oo), 

< 

u = V = 0 on r X (0, -|-oo), 

_ {u,ut){0) = {v,vt){0) = in O 

The decay speed is shown to change corresponding to the various properties of the nonlinear 
damping p{x,ut), especially, if p{x,ut) = aut {a > 0),b{x) = b > 0, the total energy 
decays exponentially. This phenomenon indicates that the velocity coupling has different 
impact compared to the displacement coupling. Being regarded as perturbation, the coupling 
through displacements is compact, while the coupling through velocities is bounded. On the 
other hand, the controllability (asymptotic stability) of coupled wave systems with general 
coefficients is closely related to the synchronization (asymptotic synchronization), according 
to the pioneer results on synchronization by Li and Rao [a US] (see also [l0|). For the above 
reasons, we focus on the question: What kind of coefficients can lead to the exponential 
stability of the general coupled wave system by velocities? 

Suppose that O is an bounded open set in with boundary T = (90. Consider the 
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following general wave system coupled by velocities 


Utt - Au + aut + /3ut = 0 

in 

Q X (0, + 00 ) 

vtt - Au + 7Ut + r]Vt = 0 

in 

Q X (0, + 00 ) 

u = V = 0 

on 

F X (0, + 00 ) 

{u,ut){0) = {v,vt){0) = 

in 

n 


where a, /3, 7 , r; e M are constants. 


( 1 . 1 ) 


Definition 1.1. System (jl.ip is said to be exponential stable if there exist constants M > 0 
and w > 0 such that 


E{t) < Me-‘^^E{0) 


( 1 . 2 ) 


where the (total) energy is defined by 

E{t) = l- f \ut\‘^+ \Vu\'^+ \vt\'^+ \Vvfdx 

^ Jn 

and uj is the corresponding decay rate. 

The task of this paper is to find the conditions that the coefficients matrix 


a /3 

.7 Vj 


(1.3) 


should 


satisfy such that System (jl.ip is exponential stable. Our main result is the following theorem: 
Theorem 1.1. System (11.11) is exponential stable if and only if 


a + ?7 > 0, 

ar] — P'j > 0 


(1.4) 


i.e., the two eigenvalues of the coefficient matrix 




both have positive real part. 


Remark 1.1. By Proposition 12.21 the equivalent form of Theorem ll.il one can see that the 
condition ()1.4p means that the two components of the solution are essentially both damped. 


The organization of the paper is as follows: In Section 2, we give the well-posedness of 
the system (tni and reduce the original general problem equivalently into the same problem 
in two canonical forms (see Proposition 12.21) . Then, we prove Proposition 12.21 in Section 3 
by three different approaches successively, that is, Riesz basis approach, multiplier method, 
frequency domain approach. Finally, some useful extension and remarks are provided in 
Section 4. 


2 Preliminaries 

For the simplicity of statements, let C? = £^(12),"Rq = and % = 

T-Lq X X Hq X . Then "H is a complex Hilbert space equipped with the inner product 


{U,U)-u= / Vu ■ Vu + yy+ 'Vv ■Vv + zz dx. 


( 2 . 1 ) 


3 




Let U = {u,y,v,z) and let Uq = {uo,ui,vo,vi). System (jl.ip can be rewritten as the 
Cauchy Problem in abstract form: 

U{0) = Uo. ( 2 . 2 ) 

where A : T>{A) —)• is defined by 


AU = (y, Art — ay — f3z, z, Av — 72 / — yz) (2.3) 

with 'D{A) = ("H^n^Q) x'Hq X ("H^n^Q) x'Hq. Obviously, A can be regarded as a bounded per¬ 
turbation to the standard wave operator, namely, A = Aq + Ai where AqU = (y. An, z, Av) 
with 'D{Ao) = T^iA) and AiU = (0, —ay — /3z, 0 , —yy — yz) with 'D(Ai) = %. 

By classical Hille-Yoshida Theorem and perturbation theory (see |20l page 76, Theorem 
1 .1]), we have 

Proposition 2.1. generates a Co-semigroup {S{t) = e*'^}t>o on 77. If Cq = (uq, ui, uq, ui) G 
77, System (jl.ip has a unique solution 7/(7) = S{t)Uo G C°([0,-|-oo); 77). If T/q G ^{A), System 
has a unique solution 7/(7) = S{t)Uo G C‘^([0,-t- 00 ); P(7l)) 0 C^([0,-|-oo); 77). Moreover, 
l|7/(7)|||^ = 2E{t) for all 7 > 0. 


For the sake of convenience of statement, we use real Schur decomposition for the coef¬ 
ficient matrix so that we can reduce the original general problem into the same problem in 
two canonical forms. 


Lemma 2.1. [9l page 79, Theorem 2.3.1] For any B = 


a /3' 

a 


G 


p2x2 


, there exists real 


orthogonal matrix P G such that BP = B where B is in one of the following two 

canonical forms: 

^ a 

-b aj 

{::) 

Moreover, ()1.4I) is equivalent to a > 0 in Case i) and to a > 0 and c > 0 in Case ii). 


i) B = 


ii) B = 


(a, 6 G M, 6 7 / 0); 
(a, 6 , c G M). 


Lemma 2.2. Let P G be a real orthogonal matrix and let {u,v) be solution of (jl.ip . 
then (u, v) = {u, v)P is solution of 


utt - Au + aut + P vt = ^ 


in fl X (0, -|-oo). 


vtt - Av + ^ ut + yvt = ^ 

u = V = 0 

^ (u,u)( 0 ) = (u°,u°)P, {ut,vt){0) = (^t^u^)P 


in 17 X (0, -|-oo), 
on F X (0, -|-oo), 
in 17 


(2.4) 


where B 



P^BP and the energy of ( 12 .4p is equivalent to that of (jl.ip . 
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Thanks to Lemmas 12.1II2.21 Theorem 11.11 is equivalent to the following Proposition 


Proposition 2.2. System (II. ip is exponential stable if and only if 
a b 


i) a > 0 when B = 

\^—5 a J 

ii) a > 0 and c > 0 when B = 


(a, 6 € M, 6 / 0); 

(a, b,c £ 


a O' 
.b c, 


3 Proof of Proposition 12.2 


3.1 Riesz basis approach 

In this subsection, we adopt Riesz basis approach to prove Proposition 12.21 in one dimen¬ 
sional case, that is = (0, tt) C M. The key ingredient is to prove the Riesz basis property 
and to analyze the spectrum of A. Additionally, we obtain the optimal and explicit decay 
rate in this situation, see Corollary 13.11 

In Case i), System dEU) can be rewritten as the Cauchy Problem ( 12.211 in one dimensional 
case, where A : 'B>{A) —)• is dehned by 

AU = {y, Uxx - ay - bz, z, v^x + by - az) ( 3 . 1 ) 

with v{A) = n nl) x x n nl) x nl. 

Let A be the eigenvalue of A and E = {u, y, v, z) G 'E>{A) be its eigenvector: 


(A/ - A)E = 0 


It is equivalent to 


Xu — y = 0 

in Til 

(A a)y -Uxx + bz = Q 

in £2 

Xv — z = 0 

in Til 

{a + X)z - Vxx - by = 0 

in £2 

Uxx + (A -|- a)Xu + bXv = 0 

in £2 

Vxx + (A -|- a)Xv — bXu = 0 

in £2 


or further 


u{0) = u{tt) = ?;( 0 ) = v{tt) = 0 
Canceling v, we obtain a forth order ODE of u: 

^xxxx 2A(A a'jUxx T [(^ T A) -1- b ]A u — 0 
u(0) = u{tt) = UxxiO) = Uxxi-^) = 0 


(3.2) 


whose general solution is given by u{x) = + A 2 e'^^^ + -(- where Vi [i = 

1, • • • ,4) are four zeroes of the algebraic equation 


— 2(A -|- a)Az/^ -|- [(a -|- A)^ -|- b^]\^ = 0 


(3.3) 
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and Ai (i = 1, • • • ,4) are constants to be determined. Then it follows that 


and consequently 


, , sinnx , . . ^. 

u{x) = -, u = ±in [n G 

n 


v{x) = + n y{x) = Au(x), z{x) = \v{x) 

— o\ 


where the eigenvalue A satisfies the characteristic equation 

A'^ + 2a\^ + {2v? + + b^)X^ + 2an^X + n'^ = 0 

or equivalently 

(A^ + (a — i 6 )A + n^)(A^ + (a + i 6 )A + n^) = 0 
Consequently, there are two classes of eigenvalues: 


^l,n ~ r,) i(=t^ „)) ^2,n ~ 9)+i(=*=^ r^) 


where 


= 


yj {a? — — 4n^)^ + 40^6^ + a? — — An? 


Y - — 

-A 77 


yielding that Xn is a decreasing function of n G It follows that 


(3.4) 

(3.5) 

(3.6) 

(3.7) 


max{5R(A±J,5R(A±J} = 5R(A+i) = 5R(A+,) = Xi - - 
nez+ ’ ’ ’ ’ 2 

Xf^ ~ Tin, A^^ ~ Tin, as n —)> +oo 

The corresponding eigenvectors are 

■ A ^tn -i . ,1 i , 

-®in = sinnx(-,-, —,-), S2n = smnx{-, -, -) 

’ n n n n ’ n n n n 

If o < 0, there exists an eigenvalue with nonnegative real part according to (I3.4p - ()3.5p . 
Then by choosing the corresponding eigenvector as the initial data, it is easy to conclude 
that System dni is unstable. If a > 0, we have 

Proposition 3.1. For any given a > 0 ,6 G M, 5 7 ^ 0, {(-F^^, forms a 

Riesz basis of the Hilbert space T-L. 

Proof. Since a > 0, 6 / 0, (11.111 has no multiple eigenvalues, thus {{EA, ^ 2 n)}n&+ 

are linearly independent. Note that 


X _ -^ _l_ 1 ^ „± _ -■ ( —, ±i, —, ::pl) (n G Z"*") 


ei„ = sinnx(-,±i,--,±l), 62 ,^ = 


sm nx I 


n n 

forms a Risez basis of % and 


n n 


4,m)'H = (e*>> = 0, Vn, m G Z+, n / m, i,j = 1,2 
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(etn> 4,n)H = {4,n^ = (ei,n> 4^^ = (ei,n> = 0> ^ 

Then using the asymptotic expansion of the eigenvalues (13.71) . we get 


4,n ~ ^tn^m + \4,n ~ -®i,nllw “ * — 1, 2 




and thus 


y~! y~! (iNi^n “ + ll®i,n “-®i,nllw) < +°0 

nSZ+ i=l,2 

By Lemma HH L^ 2 “n)}nez+ forms a Riesz basis of %. 


□ 


Thanks to Proposition 13.11 for any given Uq ^T-L, there exist 02 „)}nGZ+ 

C C'^ such that 

C0 = E E + P-8) 


ngZ+ i=l,2 


and 


E(0) = 2\\Uo\?u ~ E E (lEnl + l«*“nr) < +00 

nGZ+ j=l,2 

Moreover, the solution of (12.2p writes 


^(^)= E E Hn4’-'Eln + 4n4--'E-n) 

n€Z+ j=l,2 


(3.9) 


(3.10) 


In Case ii), consider the eigenvalue and eigenfunction: {\I — A)E = 0, € T>[A), namely 

(y, Uxx - ay, z, Vxx - by- cz) = A(u, y, v,z) £% 

Then the eigenvalue satisfies 

A"^ + (a + c)A^ + (2n^ + ac)\^ + (a + c)n^\ + = 0 

or equivalently, 


(A^ + aA + u){a + cA + n^) = 0 


thus 


aL = ': 


—a ± y/a? — AiA a 

---, if n < 


- 2’ x± 


—a ± iVdn^ — a 

- 2 -’ ^ 2 - 


Kn^< 


—c± Vc^ — 4n2 c 

-, if n < —, 

2 “2 


—c ± iV4re^ — 


(3.11) 


, if n > 


Clearly 


max 3f?(Af ) = i = 1,2 

n^Z+ ’ ’ 


Xf^ ~ Tin, A^^ ~ ±in, as n ^ +00 




'l,n 


'2,n 


(3.12) 

(3.13) 


Therefore all the eigenvalues {Af^, Xf^}n£z+ have negative real parts if and only if o > 0 
and c > 0. Certainly if a < 0 or c < 0, System ([HD is not asymptotically stable. 
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When a > 0 and c > 0, we would like to prove that the eigenvectors as well as root-vectors 
form a Riesz basis of %. For this purpose, we discuss the different situations that multiple 
eigenvalues may appear with various values of a and c. 

• Case 1. a 0 2Z+,c 0 2Z+ and a ^ c. In this case,A^^,A^^,A^^,A^^ are distinct and 
System dm) is decoupled. The eigenvectors corresponding to Af^, A^^ (re G Z+) are 

_i_ 1 b , 1 A^„ 

re re (a - c)n (a - c)n re re 

Comparing {Ef^^ with a Riesz basis of T-i: 

+ /I 

n = sin rex (-, ±1, 7-r 

re [a — c) 


zbi6 , , 1 . 

, .-)) 62^ = smrex(0,0,-,±i) 

— c)n a — c ’ re 


yields by (|3.13l) and Lemma lA.ll that {E^E^E^E^^n&+ forms a Riesz basis of 
%. Therefore, every initial data Uq G 71 can be expanded as (13.8p with ()3.9h and the 
corresponding solution is given by (|3.10l) . 

• Case 2. o G 2Z'*' or c G 2Z'*' and a ^ c. There are three sub-cases, 
o Case 2.1. a = 2m 7^ c = 2k, m,k G Z+. System (II.ip has two multiple eigenvalues: 
A^^ = —^ = —m and A^^ = — | = —k. The dimension of the eigen space corresponding to 
A^^ or A^^ is one. The corresponding eigen function corresponding to Af ^ is 

1 


Et = sinrrex (—, -1, 
’ rre 


b -b ^ 

{a — c)m 'a — c' 

Solving {A - XtmI)Er = Ef we obtain the root vector 

sinmx ^ 1 ^ b b ^ 

'^l.m o V 5 ^5 / \ 5 / 

’ 2m m (a — c)m a — c 

which is linearly independent of Ef^. Similarly , we calculate the eigenvector and root vector 


of K,k = -k- 


„4- , / 1 ^ smkx , 1 

Elk = sin kx (0,0, -1), .^2 (0,0, 1) 


By Lemma [A. 11 {Ef^,E^^,Ef^,E 2 ^}n£Z+ forms a Risez basis of Ti. For every Uq G 71 
given by (13.8|) with (|3.9p . the solution of (|2.2I) is 

Z] (3.14) 

nSZ+ i=l,2 

o Case 2.2. a = 2m,m G TA and c 0 2Z'''. A has one multiple eigenvalue = —f = 
—m. For every Uq G 77 given by (13.81) with (|3.9p . the solution of (|2.2I) is 

nGZ+ i=l,2 

o Case 2.3. c = 2k, k G lA and a 0 2Z'''. A has one multiple eigenvalue A^^ = — | = —k. 
Similar to Case 2.2, for every Uq G 77 given by (|3.8I) with (|3.9I) . the solution of (12.21) is 


!!(*)= E E + 

n£Z+ 2=1,2 


(3.16) 
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• Case 3. a = c ^ 27^^. The two classes of eigenvalues coincide Aj“„ = A^^, 

and the dimension of eigenspace varies with respect to b. 

o Case 3.1. 6 = 0. System (HI]) is decoupled. The dimension of eigenspace of Af ^ is two. 
The eigenvectors are 

1 A^ 1 A^ 

Ef = sinnx (—, ^^,0,0), Ef = sinnx (0,0, —, 

n n n n 

which forms a Riesz basis of H . For every [/q G "R given by (13.8|) with (13.Op . the solution of 
( 1 ^ is still dSTOp. 

o Case 3.2. 6/0. The dimension of eigenspase of A^^, A)|"„ is one. The eigenvector and 
root vector are 


1 


Ef = smnx ( 0 , 0 , -, 

n n 


eL = -- 


sinnx cl + 2\f^ a + 2Xf^ 1 1 

n ^ 6 A±„ ’ 6 ’ 2 A±„’~ 2 ^ 


2At 


which forms a Riesz basis of T-L. For every Uq given by (13.8p with (13.Op . the solution of 


IS 


(3.17) 

nGZ+ *=1,2 

• Case 4. a = c = 2m, m G 2Z. Then A^^ = A^^, A)|"„ = for all n G Z"*". 
o Case 4.1 6 = 0. System (|l.ip is decoupled. For n ^ m, the algebric degree of the 
eigenvalue Af „ is two. The dimension of eigenspase of Af „ is two and the eigenvectors are 

.1 Af 


Ef„ = sin nx 

're' re 


(-,—, 0 , 0 ), = sinrex( 0 , 0 ,-,^) 

re re ’ re re 


For n = m, the algebric degree of the eigenvalue A)*"^ is four. Actually A]*^^ = A 2 
and the dimension of its eigenspase is two. The eigenvectors and root eigenvectors are 


= —m. 


Efm = sin 7 rex(^,- 1 , 0 , 0 ), Ef^ = smmx{0,0, ^,- 1 ) 


1 


E7^ = 


m 

sinmx , 1 


2m ^m 


( ) 1 ) 0 , 0 ), E2 ^ — . , , 

m 2 rre rre 


( 0 , 0 ,-,!) 


For every [/q G "R given by (13.81) with (|3.9p . the solution of 


IS 




n.sZ+ *=1,2 


(3.18) 


o Case 4.2 6/0. For re / m, the eigenvalues A^^ are the same as in Case 3.2. While for 

re = m, the dimension of the eigenspace of the eigenvalue A^^ = A^^ = —m is one. We could 

solve successively the eigenvalue, hrst order, second order and the third order root vector: 

, 1 , sinrrex , 1 

Efm = sinrrex( 0 , 0 , —, - 1 ), ^ = —-( 0 , 0 , —, 1 ), 

’ rre ’ 2m m 

sinrrex ^ 1 —1 1 sinrrex ^3 —1 1 


rre^ z 6 rre’ 26 ’ 2m^ 
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For every Uq £71 given by (13.8p with ()3.9p . the solution of 


IS 


nez+ *=1,2 

n^m 




(3.19) 


2V 


3! 


Above all, we easily conclude by (13.9^ . (13.101) . (I3.14p . (f3.15p . (l3.1^ . (j3.17p . ()3.18p . (l3.19p that 

E{t) = 2\\U{t)\\'^ < ME{0)e-‘^\ Vt e [0, +oo) (3.20) 

for some positive constants M, u independent of the initial data. This is the end of Proof of 
Proposition 12.21 □ 

From the above proof of Proposition 12.21 we get the corollary concerning the decay rate. 

a b 


Corollary 3.1. i) If i? = 


—b a , 


(a > 0,6 G M, 6 / 0); then one has (13.201) with the 

(a > 0, 6 € M, c > 0), 


a O' 

,6 Cj 


optimal decay rate oo = —23?(A^;^) = —23ft(A^^). ii) If B = 
then one has 

E{t) = 2\\U{t)\\'^ < ME{0)tPe-‘^\ Vt G [0,+oo). 
with the optimal decay rate lo = —2max{$R(A^]^), 5R(A^]^)} and 


(3.21) 


p = < 


1 , if a = c = 2 and 6 = 0 ora = 27 ^c ora/2 = c; 

3, if o = c = 2 and 6/0; 

0 , else. 


(3.22) 


3.2 Multiplier method 


In this subsection, we apply the multiplier method to establish the decay estimates of the 
total energy. The key ingredient is to use an integral inequality (see Lemma [A.2p which leads 
to the exponential decay of the energy. In this subsection, the initial data are assumed to be 
all real functions. 

In Case i), System (jl.ll) is reduced to 


utt - An + aut + bvt = 0 
vtt - Av - but + avt = 0 
u = V = 0 

^ (u,ut)(0) = (u°,M^), (v,vt)(0) = (u°,u^) 


in fl X (0, +oo), 
in fl X (0, +oo), 
on P X (0, +oo), 
in Q 


(3.23) 
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Using the multiplier ut to u-equation and vt to u-equation, we obtain 

/ ut{utt - Au + aut + bvt) + vt{vtt - Av + avt - but)dx = 0 
Jn 

Integrating by parts and using the definition of the total energy E{t), we get 

dE{t) _ f Vidx < 0 

Jn 


(3.24) 


dt Jn 

which indicates that the energy decays if and only if a > 0. Then, we can follow the proof 
of |12l page 114, Theorem 8.13] to conclude by LaSalle Invariance Principle [13] that System 
Ea is asymptotically stable if a > 0. 

Using the multiplier u to u-equation and v to u-equation, we obtain for 0 < 5 < T that 

/ / u{utt — Au + aut + bvt) + v{vtt — Av — but + avt)dxdt = 0 

Js Jn 

and after integration by parts, 

T 


/ 

Jn 


uut + vvtdx 


ul + vfdxdt + 


s Jn 

T 


jVup -|- jVupdxdt 


-II 

-I-— f (u^ -|- v‘^)dx + b f [ uvt — vutdxdt = 0 

^ Jn s Js Jn 


thus 


r/ 


jVup-I-jVupdxdt = f j Ut+Vtdxdt — b f f uvt — vutdxdt 

Js Jn Js Jn 


- / UUt + VVt + 

Jn J 


dx 


Thanks to Cauchy Inequality and Poincare Inequality, we get for every e > 0 
[ f jVuP-I-jVuPdxdt < (1 -|- Ce) [ f uf + vfdxdt + e f f u^ + v^dxdt 

JsJn ~ JsJn Js Jn 

- / UUt + VVt + ^{u'^ + v‘^) 

Jn J 

< (1 + C'e) 

+ CiE{T)+E{S)) 


dx 


^ f uf + Vtdxdt + Ce f f jVup-|-jVupdxdt 

s Jn Js Jn 


Noting (I3.24P implies for a > 0 

rT 


[ [ u^t+ vldxdt = -{E{S) - E{T)) < -E{S) 

Js Jn « « 

Choosing e > 0 small we get 

[ [ jVup + jVupdxdt < 

Js Jn 
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and further 


/ 


E{s)ds < CE{S), VO < 5 < T 


(3.25) 


Thanks to Lemma lA.2[ we conclude that E{t) decays to 0 exponentially. 

In Case ii), System dEU) is reduced to 

utt — An + aut = 0 
vtt - Av + but + cvt = 0 

u = V = 0 

^ {u,ut)iO) = (u,vt)( 0 ) = 

It is easy to see that if a < 0, then the energy of u does not decay; while if c < 0, by 
taking u = 0, then the energy of v does not decay. Hence, the total energy E{t) of (I3.25|) 
does not decay to 0 if a < 0 or c < 0. It remains to prove that if a > 0 and c > 0, E{t) 
decays to 0 exponentially. 

For K > 0, we set the equivalent energy 

K|Vttp |Vup 


in X (0, +oo), 
in X (0, +oo), 
on r X (0, +oo), 
in 12 


F f , Vt 


+ 


-dx 


(3.26) 


Using the multiplier KUt to w-equation, vt to u-equation, 

dE^{t) 


dt 


= — Kaut + butVt + cvtdx 

Jn 


Choosing k > 0 suitably large, namely, 6 ^ — Akuc < 0 then there exists <5 > 0 such that 

dE^{t) 


dt 


< -6 


KUt + vfdx 


(3.27) 




Using the multiplier ku to ti-equation, vt to u-equation, 

pT I- i-T 

KUUt + VVfdx 

S 




fl 

Js Jn 

+ 


KU^ + vfdxdt + 


+ \Vv\^dxdt 
T 


- I {kqu^ + cv^)dx +b I vufdxdt = 0 
^ Jn s Js Jn 


thus, 


[ [ KlVttP + iVuPdxdt = f f Kuf + vfdxdt — b f f vufdxdt 

Is Jn JsJn Js Jn 

T 

dx 


I 

Jn 


1 


— I KUUt + VVt + -[kUU^ + CV^) 


Thanks to Cauchy Inequality and Poincare Inequality, we get for every e > 0 


f f «;|Vu|^ + \Vv\‘^dxdt 

Js Jn 


< 


f! 

's Jn 


Kuf + vfdxdt + 


fl 

Js Jn 


1 


Is Jn 

.2 I „„,2^ 


ev^ + CeKu^dxdt 


— / KUUt + VVt + -{kUU^ + CU^) 

Jn 2 

r-T 

Kuf + vidxdt 


dx 


<{l + Ce)[ [ 

Js Jn 

+ Ce [ f \S/v\'^dxdt + C{E^{T) + E^{S)) 
Js Jn 
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Similarly as as in Case i), we can prove by (), '1.271) and choosing e > 0 small the integral 
inequality for E,^: 



< CE^iS), 


yo<s <T. 


Then it follows by Lemma lA.21 E^(t) (or equivalently, E(t)) decays to 0 exponentially. This 
ends the proof of Proposition 12.21 □ 


3.3 Frequency Domain Approach 


In this subsection, we use the frequency domain approach to prove Proposition 12.21 More 
precisely, we want to get the exponential stability of the semigroup through the uniform 
estimate of the resolvent on the imaeinarv axis bv Lemma IA.3I [111 I22j Isee also m)- 

For Uq = {u^, , v^) e 71, the solution of () 2 . 2 p is U{t) = S{t)Uo where {S{t)}t>o is 

the associated Co-semigroup of operator and \\S{t)Uo\\‘^ = 2E{t) for all t >0. 

In Case (i), the energy relation (|3.24p yields that the energy decays, or equivalently 
{S{t)}t>o is a contraction, if and only if a > 0 . Next, we prove Conditions (jA.ip and (IA.2P 
in Lemma lA.31 are satisfied for a > 0. 

We start, by contradiction arguments, to assume that dATl) does not hold, i.e., there 
exists ^ G M and U = (u, y, v, z) G 'D{A) with WUW'n = 1 such that — A)U = 0, namely, 

i^tt — y = 0 in Til 

(i^ -|- a)y — Au + bz = 0 in 

< 

i^v — z = 0 in Til 

(i^ -|- a)z — Av — by = 0 in 

By definition (|2.ip . we easily calculate 


3f?((i^ - A)U, U)n = a{\\y\\l, + \\z\\%) = 0 

Then it follows that y = z = 0 in and thus 

i^u = 0 in £2 
Au = 0 in £2 

< 

i^v = 0 in £2 
Au = 0 in £2 


(3.28) 


The theory of elliptic equation with Dirichlet boundary condition implies that u = v = 0 £ 
Consequently we have U = 0 £ H which contradicts with ||t/||-^ = 1. Therefore 
(1135 is satisfied for a > 0. 

Next we continue to prove that ()A.2p holds for a > 0. Otherwise, thanks to the continuity 
of the resolvent TZ{i^, A) with respect to ^ G M, (IA.2p is not valid at oo, i.e., there exists {Un = 
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{Un,yn,Vn,Zn)}n^Z+ ^ V{A) with \\Un\\H = 1 and {U}n€Z+ C M with + 00 ,||(iCn “ 

~^)Un\\H —>• 0 as n —>■ + 00 , namely, 

i^nUn -Vn^O 111 

(i?n + a)yn - ^Un + bzn^ 0 in 

< 

^CnVn - Zn^O in "HJ 

^ (i?n + a)zn - Avn - byn ^ 0 in 

Similarly as (I3.28p . we have 

5R((i^ — A)Un, Un)'H = o(||yn||£2 + ||-2n||£2) (3.29) 

On the other hand, 


|$R((iC„ - A)Un, Un)n\ < \mn " A)Un\\H ' WUnWu = IKi^n 


•^)Un\I h 0 - 


Then —>• 0 € C^, yn ^ 0 ^ C?, and thus 


AU'fi 

i'^n 

AVn 

i'^n 


0 

0 

0 

0 


in £2 
in £2 
in £2 
in £2 


Thanks to Nirenberg Inequality [U page 135, Theorem 5.2], we get 


||VUn||£2 < CjlAu,, 


I 2 I 
\c‘^\ 


Uri 


1 

I 2 
l/:2 


<C'|| 


AUr, 


i 

I 2 
I £2 


I 


I 2 
I £2 


and in a same way Vvn —>■ 0 € Consequently ||[/n||-H —^ 0 as n —>■ +oo. This contradicts 
with ||t/n||w = 1 and implies (IA.2p for a > 0. The exponential decay of the semigroup of 
operator {S'(t)}i>o associated to ( 12 . 2 p is a consequence of Lemma [A.31 

In Case (ii), it is easy to see, as in Section 3.2, that if a < 0 or c < 0, system (13.2511 is not 
asymptotically stable. It remains to prove that (I3.25P is exponentially stable if a > 0 and 
c > 0 . 

We introduce an equivalent inner product of dZI]), 


{U,U)-}{^= / kVu ■ Vu + Kyy+ Vv ■ Vv + zz dx, (3.30) 

Jn 

where k > 0. Clearly, H is a Hilbert space under (•,-)'h„ and is denoted by Hk with the 
corresponding norm || • For Uq = (u°,u^) G "Hr, the solution of (12.211 is U{t) = 
S{t)Uo where {S{t)}t>o is the associated Co-semigroup of operator and ||5'(i)Co|||^^ = 2Ef^{t) 
for all f > 0. By Lemma lA.31 it suffices to prove that Conditions (lA.ip and (IA.2I1 hold for 
a > 0 and c > 0 . 
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We start, by contradiction arguments, to assume that does not hold, i.e., there 

exists ^ G M and U = {u, y, v, z) G V(A) with HC/H-h,; = 1 such that (i^ — A)U = 0, namely, 

i^rt — y = 0 in "Hq 

(i^ + a)y — Au = 0 in 

< 

i^v — z = 0 in "Hg 

(i^ + c)z — Av + by = 0 in 

We compute 

- A)U,U)n^ =a«:||y||^2 +c||z||^2 +b?R.{y,z)c2 =0. (3.31) 

Taking k > 0 suitably large, namely, 6^ — 4«;ac < 0 then there exists (5 > 0 such that 

0 = 5R((ie - A)U, U)h. > 5(112/11^2 + |k||^2). (3.32) 

Then y = z = 0 G C'^ and thus 

u = 0 in 

Au = 0 in 

< 

V = 0 in 

Av = 0 in 

\ 

Obviously the theory of elliptic equation with Dirichlet boundary condition implies that 
u = V = 0 G 'h? n^o- Consequently we have U = 0 G which contradicts with 11^^11-//,^ = 1 
and yields that (jA.lD is satished for a > 0 and c > 0. 

Next we continue to prove that (|A.2|1 holds for a > 0. Otherwise, thanks to the continuity 
of the resolvent Tl{i^,A) with respect to ^ G M, ()A.2D is not valid at oo, i.e., there exists 
{Un = {Un,yn,Vn,Zn)}neZ+ C V^A) with ||t4||w, = 1 and {Cn}nGZ+ C M with 
+00, ||(i^n — •^)Un\\nK —^ 0 as n —+oo, namely, 

/ 

Vn ^ 0 

“1“ y 0 

< 

bljn ^ 0 

Let K > We get similarly as (|3.31l) and (I3.32p that 

3^((*?n UAj'Hf^ = 0 >k\ Iz/n 11^2 T c| 1^2 T b'Si{yn, ^ri)Cp. > 5(| | 2 /nl l /;2 + | I'^^nl 1^2) (3.33) 

for some 5 > 0. On the other hand, 

-A)Un,Un)Hj < - A)Un\\n^ ' \\Vn\\n. 

— Il(i'^n “ •^)Un\\Hn 0- 


in "Hg 
in £2 

in "Hq 
in £2 
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Hence we have —>• 0 G ^ 0 € and thus 



^ 0 

in 

c 

AUr, 


U 

in 

L 


^ 0 

in 

£■ 

AVr, 


^ 0 

in 

L 


Thanks to Nirenberg Inequality [U page 135, Theorem 5.2], we get 

i i Au i i 

||Vu„||£2 < C||Attn||£2|kn||£2 < C'||-^||£2||CnWnll£2 0 

sn 


and in a same way ^ 0 G Consequently ||[/n||^„, —)• 0 as n —)• +oo. This contradicts 

with IIHjiII-Hk = 1 implies ()A.2p for a > 0 and c > 0. The exponential decay of the 

semigroup {S{t)}t>o associated to (|2.2I) is again a consequence of Lemma rA.31 


4 Extension and Remarks 


Theorem 4.1. If the boundary conditions in (HI]) are replaced by mixed types of boundary 
conditions, we can adopt the above three approaches to conclude that ()1.4p is still the sufficient 
and necessary condition of exponential stability of the system. For instance, 


utt - Art + aut + f3vt = 0 
vtt - Au + 7Ut + r/ut = 0 
< u = V = 0 
du dv 

, {u,ut)i0) = iv,vt)i0) = 

where Tq U Ti = dn, Tq fl Ti = 0. 


in H X (0, +oo), 
in H X (0, +oo), 
on To X (0, +oo), 
on Ti X (0, +oo), 
in H 


Theorem 4.2. If the damping and the coupling are both acted on the boundary, we can adopt 
the multiplier method to conclude that (II.4p is still the sufficient and necessary condition of 
exponential stability of the following system 


Utt - Au = 0 


in H X (0, +oo), 


Vtt - Au = 0 
u = V = 0 

du o fx 

— + aut + pvt = 0 
ou 

dv 

— + 7Ui + r]Vt = 0 

du 

{u,ut){0) = iv,vt){0) = iv°,v^] 


in H X (0, +oo), 
on To X (0, +oo), 

on Ti X (0, +oo), 
on Ti X (0, +oo), 
in H 
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where Toljri = d^l, FonTi = 0 with some suitable geometric conditions (see [6l[T2]). For 
instance, one can assume that there exists xq G M", such that Fq = {x G dQ\{x—xo)-i'{x) < 0} 
and Fi = {x G (90|(x — xq) ■ i'(x) > 0} / 0 where i'{x) is the unit outer normal hied. 

Theorem 4.3. If we consider a wave system with multiple propagation speeds 
utt — Au + aut + fivt = 0 in x (0, +oo), 

ujt - c^Au + 7Ui + ryut = 0 in nx(0,+oo), 

< 

u = u = 0 on F X (0, +oo), 

(u,ut)(0) = (u,ut)(0) = in n 

with c 7^ 1, we can adopt the multiplier method or the frequency domain approach to prove 
the exponential stability of the above system under the assumption 

a > 0, 

< r/ > 0, 

or/ — /37 > 0 

Remark 4.1 (Comparison of Three Approaches). The advantage of Riesz basis approach is 
to give the explicit expression of the solution so that the relation between the exponential 
stability of the system and the coefficients can be relatively easy to discover. In addition, the 
optimal decay rate can be obtained through very careful analysis of both the spectrum and 
the eigenvectors. However, the Riesz basis properties of the eigenvectors are hard to check in 
higher dimensional case. The multiplier method is rather simple to apply without restrains 
in space dimension and works for variable coefficients case (even for nonlinear problems), but 
it requires strong geometric assumptions. The frequency domain approach is also applicable 
for all space dimension and for the variable coefficients case without much information of the 
eigenvalues and eigenvectors. Nevertheless, the optimal decay rate can not be obtained in 
general by multiplier method or by frequency domain approach. 


A Appendix 


Lemma A.l. Let {en}nez+ be a Riesz basis on Hilbert space "H. Let {-En}nez+ be a basis 
on 7-i. If W^n — CnlP < +00, then {En\nez+ is a Riesz basis on T-L. 

n£Z+ 

Lemma [A. II is an equivalent form of Bari’s Theorem, see [SJ page 317, Theorem 2.3]. The 
special case that {en\n&z+ is an orthogonal Riesz basis of H can be found as [211 App. D, 
Theorem 3]. 


Lemma A.2. [121 page 103, Theorem 8.1] Let E : [0, +oo) —>■ [0, +oo) be a non increasing 
function Suppose that there exists C > 0 such that 

rT 

E{s)ds < CE{S), VO < 5 < T, 


/ 

Js 
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then 


E{t) <E{d)e^-i, yt>c 

Lemma A.3. [191 page 4, Theorem 1.3.2] Let {S{t) = e^-^}t>o be a Co-semigroup of contrac¬ 
tions on a Hilbert space T-L. Then {S'(t)}t>o is exponentially stable if and only the following 
conditions are satished 

(1) iM = {ie|eeM}Cp(A); (A.l) 

(2) limsup ||7^(i^, A)|| <-1-00. (A.2) 

ICK-l-oo 
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